Resistivity of Mixed-Phase Manganites 
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The resistivity pdc of manganites is studied using a random-resistor-network, based on phase- 
separation between metallic and insulating domains. When percolation occurs, both as chemical 
composition and temperature vary, results in good agreement with experiments are obtained. Sim- 
ilar conclusions are reached using quantum calculations and microscopic considerations. Above the 
Curie temperature, it is argued that ferromagnetic clusters should exist in Mn-oxides. Small mag- 
netic fields induce large pdc changes and a bad-metal state with (disconnected) insulating domains. 

PACS numbers; 71.10.-w, 75.10.-b, 75.30.Kz 
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The study of manganites is one of the main areas of re- 
search in Strongly Correlated Electrons ||l|. Three main 
reasons have triggered this wide interest: (1) The low- 
bandwidth materials have unexplained transport proper- 
ties. They are insulators at room temperature, chang- 
ing into bad metals at low temperatures. A sharp peak 
in the resistivity pdc appears at the ferromagnetic (FM) 
transition. Small magnetic fields turn the insulator into a 
metal, with a "Colossal" Magneto- Resistance (CMR). (2) 
The phase diagram T-x (T=temperature, x=hole den- 
sity) is rich, with complex spin, charge, and orbital order. 
(3) Mn-oxides have intrinsic inhomogeneities in most of 
the T-x plane even in single crystals ||,^. 

This challenging behavior has been addressed by pre- 
vious theoretical studies. Regarding item (2), the var- 
ious FM, antiferromagnetic (AF), orbital-ordered, and 
charge-ordered (CO) phases have already emerged from 
simulations and mean-field approximations y. Regard- 
ing item (3), phase separation (PS) has been proposed to 
explain the inhomogeneities. PS can be: (a) electronic, 
with nanometer-size clusters ^, or (b) structural, where 
disorder can induce up to micrometer-size clusters and 
percolation, when influencing on first-order transitions 
Q. However, the explanation of transport, item (1), is 
more complicated since estimations of pdc are notoriously 
difficult. In addition, the prominent inhomogeneities of 
Mn-oxides [|| have not been incorporated into realistic 
Pdc calculations p]. The behavior of pdc remains unex- 
plained, although it is central to manganite physics. 

Our goal in this paper is to present a rationalization 
of the Pdc vs T curves of Mn-oxides based on the cur- 
rently prevailing phase-separated/percolative framework 
for these compounds. In this context, items (1) and 
(3) above are closely related. The analysis necessarily 
involves phenomenological considerations, since percola- 
tion cannot be addressed on sufficiently large lattices us- 
ing accurate microscopic models. However, the /x-meter 
clusters in experiments H strongly suggest that a coarse- 
grain approach should be sufficient. In addition, results 
of microscopic calculations presented below are consis- 



tent with those of the macroscopic approach. 

The main concept introduced here is summarized in 
Figs.la-b. The manganite state in the CMR regime is 
assumed to be percolative, with metallic filaments across 
the sample (Fig. la). Percolation indeed occurs in mod- 
els g, and in many experiments |Q. The insulating and 
metallic (percolative) regions are assumed to have resis- 
tances Ri(T) and R^'^(T), respectively, as sketched in 



Fig. lb. R^ is large at T=0 due to the complex shape 
of the conducting paths, and grows with T as in any 
metal, eventually diverging when the percolative path 
melts with increasing T. Note that at room temperature 
Ri<R^' and, thus, most of the conduction in this regime 
occurs through the insulator. On the other hand, Ri is 
so large at low-T that current can only flow through the 
percolative paths. This suggests a simple two parallel 
resistances description (Fig. lb), where a peak in the ef- 
fective resistance at intermediate T is natural. 

To substantiate this idea, flrst consider results ob- 
tained using a random-resistor-network that mimics the 
prominent FM-CO mixtures [|| in Mn-oxides in the CMR 
regime. Two-dimensional (2D) and three-dimensional 
(3D) square and cubic clusters are used, with link re- 
sistivities randomly selected as metallic (pm) or insulat- 
ing (p/), with a flxed metallic fraction p (which in, e.g., 
(La5/8„yPry)Ca3/8Mn03 (LPCMO) it is proportional to 
the amount of La [||). The lattice spacing of this effec- 
tive network is comparable to the FM or CO domain size, 
much larger than the Mn-Mn distance. The actual val- 
ues of Pm and pi vs T were directly taken from LPCMO 
data 1^ (y=0.00 and 0.42, respectively), and for simphc- 
ity they are used in both 3D and 2D clusters. Other 
materials were tried and the analysis below does not de- 
pend qualitatively on the reference compounds. 

The Kirchoff equations for the network were solved it- 
eratively on large clusters using well-known techniques 
0, and the net resistivity was found. Typical results are 
shown in Fig.lc. Only the limiting cases p=0 (all insu- 
lator) and 1 (all metal) are taken from experiments. As 
expected, a percolative regime exists between p=0.4 and 



0.5, where pdc{T—Q) is as large as in LPCMO and other 
materials ^. pdc has insulating behavior at room-T, even 
for p as high as 0.65, while at low-T a (bad) metallic be- 
havior is observed. A broad peak appears at intermediate 
T's and p's. Similar results exist in 3D (inset of Fig.lc). 
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FIG. 1. (a) Mixed-phase state near percolation. Arrows 
indicate conduction through insulating or metallic regions 
depending on T (see text), (b) Two-resistances model for 
Mn-oxides. Effective resistance Rcff vs T (schematic) arising 
from the parallel connection of metallic (percolative) R^"^ s-nd 
insulating Ri resistances, (c) Net resistivity pdc of a 100x100 
cluster vs T, at the indicated metallic fractions p. Inset: Re- 
sults for a 20^ cluster with (from top) p=0.0, 0.25, 0.30, 0.40 
and 0.50. In both cases, averages over 40 resistance config- 
urations were made (errors the size of the points). The p—1 
and limits are from Ref. ^. Results on 200x200 clusters 
(not shown) indicate that size effects are negligible. 

It is remarkable that Fig.lc is already in good quali- 
tative agreement with some Mn-oxide experiments, such 
as for Lao.96-yNdyKo.o4Mn03 Q, or even non-manganite 
materials, such as CaFci-xCoxOa where an AF-FM com- 
petition occurs iQ. However, many manganites present 
a more pronounced pdc peak at intermediate T's. To 
reproduce this feature, it is necessary to introduce a per- 
colative process not only as p (or x) varies, but also as 
T changes. This is reasonable since the metallic compo- 
nent triggered by ferromagnetism is sensitive to T, and 
the FM clusters shrink in size as T increases. This pro- 
posal was tested qualitatively using two models: (i) the 
Random Field Ising model (RFIM), that describes the 
disorder-induced PS ||^, with spin up and down crudely 
representing the competing metal and insulator, and (ii) 
the well-known one-orbital model (with parameters t, 
Jh, and J' representing the Cg-hopping, Hund coupling, 
and Heisenberg exchange among t2g-spins, respectively). 
The latter is supplemented by a term J2i4>i'^ii with (j)i 
randomly taken from [-W,W] and Uj the on-site density 
at site i. This disorder generates coexisting clusters near 



first-order FM-AF transitions, as explained in Ref. |^. 

In Fig. 2a, it is shown a portion of a typical Monte Carlo 
(MC) simulation of the RFIM at T=1.6 on a 500x500 
cluster, with fixed random fields taken from [-1.0,1.0] 
(J=:l is the FM Ising coupling). These parameters are 
the same as in Ref. H. Three main domains were found: 
spins- up (black), spins-down (white), and regions with a 
small spin expectation value (grey) . The generation with 
T of these paramagnetic (FM) areas in the surface of the 
up and down domains weakens the percolative tendencies 
of the RFIM: for example, in Fig. 2a domains connected 
at T=0, become disconnected at finite-T. Similar behav- 
ior occurs in the microscopic one-orbital case, as shown 
in Fig. 2b for a disordered configuration with percolative 
features lO]. Increasing T from ~0.0 to 0.05t decouples 
the two FM regions. FM, AF, and PM regimes dominate 
at finite T, as in the RFIM. The T=0.05t Drude weight 
(not shown) is much smaller than at T'^O.O. 
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FIG. 2. (a) MC results for the Random Field Ising Model 
at T=1.6 (J=l), and a fixed set of random fields taken from 
[-1.0,1.0]. Shown is a 100x100 subset of a 500x500 cluster 
with periodic boundary conditions (PBC). Black, white, and 
grey denote regions with Ising spin mean-value (a) in the in- 
tervals (0.66,1.0), (-1.0,-0.66), and (-0.66,0.66), respectively. 
The arrows indicate regions where percolation occurs as T de- 
creases, (b) MC results for the one-orbital model on a PBC 
8x8 cluster, Jh=8, J'=0.0, t=l, density {n)=0.875, and a 
(fixed) realization of disorder (W=0.5) with percolative char- 
acteristics at T~0. Black and white denote FM and AF re- 
gions, respectively. Grey domains are FM at T~0, but the 
nearest-neighbor (NN) spin correlations become much smaller 
(PM) at T=0.05. (c) Net pdc of a 100 x 100 cluster as in Fig.lc, 
but with p changing with T (representative values indicated) . 
Averages over 40 resistance configurations are shown. 

To incorporate the indications of T-induced percola- 
tion (Figs.2a-b) in the phenomcnological approach, a T- 
dependent metallic fraction p(T) is needed, p must de- 
crease as T grows, should vary rapidly near the Curie 



temperature Tq as the niagnetization does, but other- 
wise its T-dependence is unknown. Fortunately, the qual- 
itative results using several functions are similar, and a 
typical case is shown in Fig. 2c. The pdc^ obtained by 
this simple procedure now clearly resemble those found 
in experiments, with a robust peak at intermediate T's, 
moving toward lower T's as the system becomes more in- 
sulating. This agreement with experiments is unlikely to 
be accidental, and justifies a posteriori our assumptions. 
Note that if our approach is correct, consistency requires 
that in Mn-oxides above Tc there should exist (discon- 
nected) FM clusters on an insulating matrix, since p does 
not drop abruptly to at Tg. A new temperature scale 
T* is predicted, with those FM clusters existing in the 
range Tc<T<T*. The density of states in this regime 
likely has a pseudogap,_ba,sed on previous investigations 
of mixed-phase states |^^ . 

Consider now nanometer-scale clusters. Here quantum 
effects cannot be neglected. However, the problem is still 
too difficult to be treated microscopically, and an effec- 
tive description is needed. For this purpose, instead of 
a resistor network, a 3D lattice model with NN electron 
hopping (and zero chemical potential) is here used, with 
link hopping amplitudes randomly selected to be either 
"metallic" (tu) or "insulating" (tj), representing effec- 
tive hoppings through the nanoclusters with 1 (0) cor- 
responding to the FM (AF) regions of the microscopic 
model at large Jjj. Such "lattice of quantum wires" has 
been used before to study quantum percolation |I2] . The 
cluster conductance C (in e^/h units) is calculated using 
the Kubo formula within a Landauer setup |l3] , connect- 
ing the finite clusters to semi-infinite leads using an in- 
finitesimal voltage drop. The self-energy of these leads is 
known exactly |l3], and the lead on the, e.g., right is used 
for an iterative calculation of the self-energy from right 
to left through the cluster, supplemented by a self-energy 
matching at the left end. With the Green function ob- 
tained by this procedure, C is evaluated using known for- 
mulas p^ . Since this is quantum percolation, the study 
is only in 3D (2D localization leads to zero conductivity). 

The hoppings tjvi and ti are not available from exper- 
iments. However, tM should decrease with increasing T 
following the FM-phase magnetization, while ti increases 
with T (since, e.g., the zero-conductivity T=0 AF con- 
figuration disorders as T grows). C was obtained on up 
to 20^^ clusters using the tM(T) and ti(T) in the inset of 
Fig. 3a, but the results do not depend qualitatively on the 
particular functions used, as long as tj changes rapidly 
with T near room-T, as pi does in experiments. The 
metallic fraction p in Fig. 3a was made T-dependent as 
in Fig. 2c, and the critical percolation at T=0 is expected 
to be located near pc~0.45 jl^]. Results are shown in 
Fig. 3a [|^. Once again a remarkable qualitative agree- 
ment with experiments is obtained, suggesting that both 
nano- and micro-meter clusters lead to similar results. 
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FIG. 3. (a) Inverse conductance of a tight-binding (NN 
sites) 20^ (L=20) cluster with hoppings tM and tj. Tempera- 
ture is considered through changes in those hoppings, shown 
in the inset, p is the fraction of metallic links, and represen- 
tative values are indicated. The chemical potential was set 
to 0. (b) Inverse conductance of the half-doped one-orbital 
model on a PBC 8x8 (L=8) cluster, and Jh^co (t=l), vary- 
ing J'. Disorder strength is A=0.1. Averages over 15 disorder 
configurations are shown. Error bars in Figs. 3b and 4a are of 
the order of the oscillations of each curve. 

For completeness, C was also calculated using micro- 
scopic models on small lattices. A MC simulation of 
the one-orbital model on an 8x8 PBC cluster and den- 
sity x=0.5 was performed. From previous work g, it 
is known that a metal-insulator first-order transition oc- 
curs at J^~0.07 (if Jh=cx)). Disorder is introduced such 
that in the NN-sites link (ij) the hopping is tij=l-|-(5ij 
and Heisenberg coupling is J(j=J'(l-|-(5y), where 5^ is 
randomly taken from [-A,A], and J' is uniform. This dis- 
order makes the transition continuous 0. The MC pro- 
cedure generates t2g-spin configurations from which NN- 
sites effective hoppings can be calculated (as in double 
exchange models). These hoppings are used to evaluate 
C (Fig. 3b). Due to the disorder, C interpolates smoothly 
from metal to insulator varying J' (otherwise a discon- 
tinuous transition occurs). C^^(T~0) can be very large, 
but finite, if the appropriate value of J' is selected [|5|. 
In this respect the result has clear similarities with those 
of the macroscopic approach. However, the full shape of 
the experimental pdc curves is difficult to reproduce with 
microscopic models on small clusters where percolation 
cannot be studied. Nevertheless, for the (few) disorder 
configurations with percolative-like characteristics found 
on small systems (as in Fig. 2b), the associated C vs T 
has a broad maximum at a finite T. 

Consider now nonzero magnetic fields (/i). In the 
random-network model mimicking coexisting FM-AF re- 
gions, a small h will increase the FM fraction p by a 
concomitant small amount. However, near percolation 
tiny modifications in p can induce large pdc changes, as 



shown in Fig. 2b where a 5% modification in p at low-T 
can alter pdc by two orders of magnitude. In the percola- 
tive regime, "small" perturbations can drastically change 
the conductivity. This analysis predicts that the metallic 
state reached from the insulator with magnetic fields is 
not homogeneous but must still have a substantial frac- 
tion of insulating clusters. This is consistent with the 
experimental large Pdc(T=0) of such a state. 
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FIG. 4. (a) Inverse conductivity of the half-doped 
one-orbital model on a 64-site chain, with Jh=cxd, J'=0.14, 
t=l, and A=0.03, varying a magnetic field as indicated. The 
data shown corresponds to one disorder configuration, but 
results with other configurations are similar, (b) Effective 
resistivity of a 100x100 network of resistances. Results at 
A(T=0.0 (full circles, open triangles, and open squares start- 
ing at T=0 with p=0.45, 0.50 and 0.70, respectively) are the 
same as found in Fig. 2c. Full triangles, inverse open triangles, 
and diamonds, correspond to the same metallic fractions, but 
with a small addition to the insulating conductivity {Aa=0.1 
(ficm)^^), to simulate the effect of magnetic fields (see text). 

Another effect contributes to this phenomenon. It ex- 
ists even on chains where percolation does not occur, 
and it is illustrated in Fig. 4a where C~^ is shown using 
the microscopic half-doped one-orbital model at J'=0.14 
where the system is at a FM(metal)-AF(insulator) tran- 
sition, the latter with the periodic spin structure ttiii 
as shown in Ref. M. The field is introduced as hj^i Mf, 
where Mf =sf + (3/2)Sf , with Sf the z-component of the 
classical t2g-spin at site i with norm 1, and sf the spin 
of the mobile electron at the same site. Disorder in the 
hopping and J' (as in Fig. 3b) of strength A=0.03 pro- 
duces coexisting FM-AF clusters B. The nearly perfect 
AF-links at low-T induce a huge C^^ at h—Q. However, 
field modifications of just O.Olt (-17 Tesla, if t=0.2eV) 
produce dramatic changes in C~^ at low-T (Fig. 4a). 
The resulting C~^ at h^O.O is still large, but compared 
with C~^{h—0) the effect is notorious. An analysis of 
the spin correlations vs h shows that these large resis- 
tance changes mainly originate in the AF regions, since 
small fields produce a small spin canting and concomi- 



tant small conductivity, creating a valve effect between 
metallic domains. In real manganites, a relatively mod- 
est Pdc change in the insulating regions could contribute 
appreciably to the large MR. To simulate this effect, 
Pi of the 2D network of Fig. 2c was slightly modified 
{Aadc=0-^{^cm)~^) , with pM untouched. The resulting 
Pdc changes (Fig. 4b) are indeed large at low-T, compara- 
ble to those obtained changing p by a few percent. 

Summarizing, the manganite pdc was studied within 
the PS-framework using a semi-phenomenological ap- 
proach. At room-T, conduction predominantly through 
the insulating regions leads to dpdc/dT<0, while at low-T 
the metallic filaments carry the current. The large MR 
induced by magnetic fields is caused by small changes 
in the metallic fraction p and/or in the conductivity of 
the insulator, effects which severely affect transport near 
percolation. Our approach provides a simple explanation 
of the CMR effect, without invoking polaronic or Ander- 
son localization concepts, and independently of the origin 
(Coulomb vs Jahn- Teller) of the competing phases. 
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